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The D~ meson can decay purely leptonically via an- 
nihilation of the c and s quarks into a W~ boson [l|. 
In the Standard Model (SM), the leptonic partial width 
T(Dj— >£~Di) is given by 
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where Md s and rri£ are the D~ and lepton masses, re- 
spectively, Gp is the Fermi coupling constant, and V cs 
is an element of the Cabibbo-Kobayashi-Maskawa quark 
mixing matrix. These decays provide a clean probe of 
the pseudoscalar meson decay constant fo 3 ■ 

Within the SM, /£> s has been predicted using sev- 
eral methods [2|; the most precise value by Follana 
et al. uses unquenched LQCD calculations and gives 
/d s = (241±3) MeV. Currently, the experimental values 
are significantly larger than this theoretical prediction. 
The Heavy Flavor Averaging Group combines the CLEO- 
c, Belle and BABAR measurements and reports fo,= 
(254.6 ± 5.9) MeV % Models of new physics (NP), in- 
cluding a two-Higgs doublet [4] and leptoquarks [5J , may 
explain this difference. In addition, fo 3 measurements 
provide a cross-check of QCD calculations which predict 
the impact of NP on B and B s meson decay rates and 
mixing. High precision determinations of fu 3 , both from 
experiment and theory, are necessary in order to discover 
or constrain effects of NP. 

We present absolute measurements of the branching 
fractions of leptonic Dj decays with a method simi- 
lar to the one used by the Belle Collaboration 0, 
An inclusive sample of D~'s is obtained by recon- 
structing the rest of the event in reactions of the kind 
e+er ^cc^DKXD*-, where D*~ -> D77. Here, D 
represents a charmed hadron (D°, D + , D*, or A£), K 
represents the K® or K + required to balance strangeness 
in the event, and X represents additional pions produced 
in the cc fragmentation process. When the charmed 
hadron is a an additional anti-proton is required to 
assure baryon number conservation. No requirements 
are placed on the decay products of the D~ so that 
the selected events correspond to an inclusive sample. 
The 4- momentum of each D~ candidate, p r , is measured 
as the difference between the momenta of the colliding 
beam particles and the fully reconstructed DKX'j sys- 
tem: p r = p e + + p e - - p D - p K - p x - p 1 - 
The inclusive D~ yield is obtained from a binned fit to 
the distribution in the recoil mass m r (DKXj) = \f$.. 
Within this inclusive sample, we determine the fraction 
of events corresponding to Dj — >/!~i> M , D~ — >e~P e , and 
D~ -^t~D t decays. In the SM, ratios of the branching 
fractions for these decays are e~ v e :^~ v^:t~ v t =2 x 10~ 5 : 
1 : 10, due to helicity and phase-space suppression. 

The analysis is based on a data sample of 521 fb _1 , 
which corresponds to about 677 million e + e~ — > cc 
events, recorded near y/s = 10.58 GeV by the BABAR 
detector at the SLAC PEP-II asymmetric-energy col- 
lider. The detector is described in detail in Refs. @, Q. 
Charged-particle momenta are measured with a 5 layer, 
double-sided silicon vertex tracker (SVT) and a 40 layer 
drift chamber (DCH) inside a 1.5 T superconducting 
solenoidal magnet. A calorimeter consisting of 6580 
CsI(Tl) crystals (EMC) is used to measure electro- 
magnetic energy. Measurements from a ring-imaging 
Cherenkov radiation detector, and of specific ionization 
(dE/dx) in the SVT and DCH, provide particle identi- 



fication (PID) of charged hadrons. Muons are mainly 
identified by the instrumented magnetic flux return, and 
electrons are identified using EMC and dE/dx informa- 
tion. The analysis uses Monte Carlo (MC) events gen- 
erated with EvtGen and JETSET [lfl, and passed 
through a detailed GEANT4 [l2| simulation of the detec- 
tor response. Final state radiation from charged particles 
is modeled by PHOTOS QJ]. Samples of MC events for 
e + e _ annihilation to qq (q = u, d, s, c, b) (generic MC) are 
used to develop methods to separate signal events from 
backgrounds. In addition, we use dedicated samples for 
D~ production and leptonic decays (signal MC) to de- 
termine reconstruction efficiencies and the distributions 
needed for the extraction of the signal decays. 

We reconstruct D candidates using the following 
15 modes: D° -> A^7r+(7r°), K-tt + tt-tt+(tt°), or 
A07r+7r-(7r°); D+ -> K-tt+tt+ (tt°), A°7r+(7r°), or 
A°7r + 7r-7r+; and A+ -> pK-TT + (n ), pK° s , or pK Q s n~TT+. 
All 7r°'s and A°'s used in this analysis are reconstructed 
from two photons or two oppositely charged pions, re- 
spectively, and are kinematically constrained to their 
nominal mass values [l4| • The in a D candidate must 
have a flight distance from the e + e~ interaction point 
(IP) greater than 10 times its uncertainty. For each D 
candidate we fit the tracks to a common vertex, and for 
each mode, we determine the mean and a of the recon- 
structed signal mass distribution from a fit to data. We 
then simultaneously optimize a set of selection criteria to 
maximize Sj V S + B, where S refers to the number of D 
candidates after subtraction of the background B within 
a mass window defined about the signal peak. Where B 
is estimated from the sideband regions of the mass dis- 
tribution. In addition to the size of the mass window, 
several other properties of the D candidate are used in 
the optimization: the center-of-mass (CM) momentum 
of the D, PID requirements on the tracks, the probabil- 
ity of the D vertex fit, and the minimum lab energy of 
7T° photons. The CM momentum must be at least 2.35 
GeV/c in order to remove B meson backgrounds. After 
the optimization the relative contributions to the total 
signal sample are 74.0% D°, 22.6% D+, and 3.4% A+ . 
Multiple candidates per event are accepted. 

To identify D mesons originating from D* decays we 
reconstruct the following decays: D* + -» D°ir + , D*° — > 
L>°7r°, D*+ -> D+tt° : and D*° -> D°j. The photon 
energy in the laboratory frame is required to exceed 30 
MeV for vr -> 77 and 250 MeV for L>*° -> L>°7 decays. 
The 77 invariant mass must be within 3 sigma of the 
7T° peak. For all D* decays, the mass difference m(D*) — 
m(D) is required to be within 2.5 sigma of the peak value. 

A K candidate is selected from tracks not overlapping 
with the D candidate. PID requirements are applied to 
each K + candidate, and a A° candidate must have a 
flight distance greater than 5 times its uncertainty. 

An X candidate is reconstructed from the remaining 
7r ± 's and 7r°s not overlapping with the DK candidate. 
In the laboratory frame, a tt^ must have a momentum 
greater than 100 MeV/c and each photon from a 7r° decay 



5 



must have energy greater than 100 MeV. We reconstruct 
X modes without 7r°'s with up to three charged pions, and 
modes with one 7r° with up to two charged pions. The 
total charge of the X candidate is not checked at this 
stage. 

Finally, we select a 7 candidate for the signal D*~ de- 
cay by requiring a minimum energy of 120 MeV in the 
laboratory frame, and an angle with respect to the di- 
rection of the D candidate momentum in the CM frame 
greater than 90 degrees. This photon cannot form a ir° or 
77 candidate when combined with any other photon in the 
event. In addition, the cluster must pass tight require- 
ments on the shower shape in the EMC and a separation 
of at least 15 cm from the impact of any charged particle 
or the position of any other energy cluster in the EMC. 

Only DKX-f candidates with a total charge of +1 are 
selected to form a right-sign (RS) sample, from which 
we extract the D~ signal yield. The charm and strange 
quark content of the DKX must be consistent with re- 
coiling from a D~ . The RS sample includes candidates 
for which consistency cannot be determined due to the 
presence of a K®. We define a wrong-sign (WS) sample 
with the same charge requirement above, but by requiring 
that the charm and strange quark content of the DKX 
be consistent with a recoil from a . The WS sam- 
ple contains a small fraction of signal events due mainly 
to DKX candidates for which the total charge is misre- 
constructed. The generic MC shows that the WS sam- 
ple, after subtraction of the signal contribution, correctly 
models the backgrounds in the RS sample. 

A kinematic fit to each DKX candidate is per- 
formed in which the particles are required to originate 
from a common point inside the IP region, and the D 
mass is constrained to the nominal value [14]]. The 4- 
momentum of the signal D*~ is extracted as the missing 
4-momentum in the event. We require that the D*~ can- 
didate mass be within 2.5a of the signal peak. For MC 
signal events, the mean is found to be consistent with the 
nominal value and a varies between 37 and 64 MeV/c 2 
depending on the number of pions in X. 

We perform a similar kinematic fit with the signal 7 
included and with the mass recoiling against the DKX 
constrained to the nominal D*~ mass 14| in order to 
determine the D~ 4-momentum. We require that the 
D~ CM momentum exceed 3.0 GeV/c, and that its mass 
be greater than 1.82 GeV/c 2 . After the final selections, 
there remain on average 1.7 D~ candidates per event, 
due mainly to multiple photons that can be associated 
with the D*~ decay. In order to properly count events 
in the fits described below, we assign weight 1/n to each 
Dj candidate, where n is the number of Dj candidates 
in the event. 

We define n x and n x to be the number of recon- 
structed and true pions in the X system, respectively. 
The efficiency for reconstructing signal events depends 
on n x . However, the n x distribution is expected to dif- 
fer from the MC simulation due to inaccurate fragmen- 
tation functions used by JETSET. To correct for these 



inaccuracies, we extract the Dj signal yields from a fit 
to the two-dimensional histogram of m r (DKX"f) versus 
n x . The PDF for the signal distribution is written as a 
weighted sum of the MC distributions for j = n x , 



S(m,n H x ) 



(2) 



3=0 



The weights Wj have to be extracted from this fit. To 
constrain the shape of the weights distribution, we intro- 
duce the parameterization Wj cx (j — a)^ e~ 7J together 
with the condition ' s p 4 . it?,- = 1. This parametrization is 
motivated by the distribution of weights in the MC. The 
value a = —1.32 is taken from a fit to MC, whereas /3 
and 7 are determined from the fit to data. 

The RS and WS samples are fitted simultaneously to 
determine the background. The fit to the WS sample 
uses a signal component similar to that used in the RS 
fit, except that due to the small signal component, the 
weights are fixed to the MC values and the signal yield is 
determined from signal MC to be 11.8% of the RS signal 
yield. The shapes remaining after the signal component 



is removed from the WS sample, Bi(m) (i 



•x 



), are 



used to model the RS backgrounds. A shape correction 
is applied to Bq to account for a difference observed in the 
MC. We add these components with free coefficients (&;) 
to construct the total RS background shape: B(m, n x ) = 
Si=o biBi(m)5(i—n x ). Thus in addition to j3, 7, and the 
total signal yield, there are 3 additional free parameters 
h(i = 0, 1,2) in the RS fit. 

Figure [T] shows the data and the results of the fit, and 
Fig. [2] shows the total RS and WS samples. The fit finds 
a mininum x 2 /ndf — 216/182 and the fitted parameter 
values are p = 0.27 ± 0.17 and 7 = 0.28 ± 0.07. These 
are different from the MC values /3 = 3.38 and 7 = 1.15 
since there are more events at low values of n x than in 
the MC. 

Having constructed the inclusive Dj sample, we pro- 
ceed to the selection of D~ — > events within that 
sample. We use the m r (DKXj) range between 1.934 
and 2.012 GeV/c 2 , which contains an inclusive Dj yield 
(N Ds ) of (67.2 ± 1.5) x 10 3 . We require that there be ex- 
actly one more charged particle in the remainder of the 
event, and that it be identified as a /i~ . In addition, we 
require that the extra neutral energy in the event, E extra , 
be less than 1.0 GeV; E extra is defined as the total energy 
of EMC clusters with individual energy greater than 30 
MeV and not overlapping with the DKX-f candidate. 
Since the only missing particle in the event should be the 
neutrino we expect the distribution of E extra to peak at 
zero for signal events. We determine the 4-momentum 
of the 9^ candidate through a kinematic fit similar to 
that described earlier in the determination of the D~ 
4-momentum, but with the /i~ included in the recoil sys- 
tem. In this fit we constrain the mass recoiling against 
the DKX-f system to the nominal value for the D~ 
To extract the signal yield, we perform a binned 
maximum likelihood fit to the m^.(DKX^[i) distribu- 
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FIG. 1: (color online) m r (DKX^f) distributions for each 
n x value. The points are the data. The open histogram is 
from the fit described in the text. The solid histogram is the 
background component from the fit. The vertical lines define 
the region used in the l~vt selections. 
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FIG. 2: (color online) m r (DKXy) distribution for the total 
WS (left) and RS (right) samples. 
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(color online) Fitted distributions 



of (a) 

m^(DKXjfi), (b) ml(DKXje), (c) E extra for D~ -^t~ uv v t , 
(d) E extra for D~ -+Tp Vv v T candidates, and (e) m(KKir'y). 
In each figure, the points represent the data with statistical 
error bars, the open histogram is from the fit described in the 
text, and the solid histogram is the background component 
from the fit. 



tion using a signal PDF determined from reconstructed 
signal MC events that contain the signal decay chain 
D*- D-j with D- yTVy,. The background PDF 
is determined from the reconstructed generic MC events 
with signal events removed. The fit is shown in Fig. |3fa) , 
and the number of signal events extracted, N^, is listed 
in Table [J 

The D~ —>fj,~Pn branching fraction is obtained from: 



N, 



N, 



fJLU 



N, 



(3) 



where the D~ —> [I'v^ reconstruction efficiency, is 
determined using the signal MC sample with j = n^, and 
e 3 D is the corresponding inclusive Dj reconstruction ef- 
ficiency. The efficiency ratios el/e^ decrease from 87% 
to 33% as j increases from to 6. The weighted aver- 
age, e^j,, and the value determined for B(Dj — > /U - ^) 
are listed in Table |U The statistical uncertainty includes 
contributions from Njj a , , and (with correlations 
taken into accounted) . The systematic uncertainty is de- 
termined by varying the parameter values in the inclu- 
sive D~ fit which were fixed to MC values, by varying 



the resolution on the Dj signal PDF (for both mass and 
n^), and by estimating how well the MC models the 
non-peaking component of the signal PDF observed in 
Figs. [T] and [2j The non-peaking signal component in the 
m r (DKX^j) distribution arises from DKX'y candidates 
in events that contain the signal decay D* ~ — > Dj'y, but 
for which the photon candidate is mis-identified and is 
due to other sources such as tt° or r\ decays, or tracks 
or K® interacting in the calorimeter. Uncertainties are 
assigned for possible mismodeling of the signal or back- 
ground rrir(DKX-f[i) distributions due to possible dif- 
ferences in the position or resolution of the mass distri- 
bution, or mismodclings of different D~ decays. Uncer- 
tainties in the efficiencies due to tracking and /i _ identi- 
fication are included. This measurement supersedes our 
previous result [l5j . 

Using a procedure similar to that for D~ — > jJTVp, 
we search for D~ — > e~v e events. The fit to the 
mf(DKX'je) distribution, shown in Fig. [3Jb), gives a 
signal yield -/V el/ consistent with 0. We obtain an upper 
limit on B(Dj — > e~v e ) by integrating a likelihood func- 
tion from to the value of B(D~ — > e~D e ) correspond- 
ing to 90% of the integral from to infinity. The like- 
lihood function consists of a Gaussian function written 
in terms of the variable BNr>^e ev with mean and sigma 
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TABLE I: Average efficiency ratios, signal yields, branching fractions, and decay constants for the leptonic D s decays. The 
first uncertainty is statistical and the second is systematic. 



Decay 






e 


Signal Yield 




fr>. (MeV) 








70.5% 


6.1 ± 2.2 ± 5.2 


< 2.3 x 10~ 4 at 90% C.L. 










67.7% 


275 ± 17 


(6.02 ± 0.38 ± 0.34) xl0~ 3 


265.7 ± 8.4 ± 7.7 


D-^t- 


~v r (r 


-> e Uef-r) 


61.6% 


408 ± 42 


(5.07 ± 0.52 ± 0.68)xl0~ 2 


247 ± 13 ± 17 








59.5% 


340 ± 32 


(4.91 ± 0.47 ± 0.54) xlO" 2 


243 ± 12 ± 14 



set to N ev and its total uncertainty, respectively. To 
account for the uncertainties on ND a s ev , the main Gaus- 
sian is convolved with another Gaussian function cen- 
tered at the measured value of Nd b £ e v with sigma set to 
the Njj a e el , total uncertainty. The value obtained for the 
upper limit is listed in Table HI 

We find D~ — > t~v t decays within the sample of in- 
clusively reconstructed D~ events by requiring exactly 
one more track identified as an e~ or from the 
decay t~ — > erv e v T or r~ — > [i~v tL v T . We remove 



struction was not applied) the background was found to 
be linear in m(KKir"/). From a fit to the m(KKir"f) dis- 
tribution, shown in Fig.[3je), we determine a signal yield 
of N KK k = 1866 ± 40 events. 

We compute the D~ — > K~K + ir~ branching frac- 
tion using Eq. ((3]). The efficiency for reconstructing sig- 
nal events is determined from the signal MC in three 
regions of the K~K + tt~ Dalitz plot, corresponding to 
<f>ir~ , K~K*°, and the rest. A variation of ^8% is ob- 
served across the Dalitz plot, leading to a correction 



events associated with D s — > [i decays by requiring factor of 1.016 on e 
mj(DKX-fn) >0.5 GeV 2 /c 4 . Since Dj ^t~v t events 
contain more than one neutrino we use E extra to ex- 
tract the yield of signal events; these are expected to 
peak towards zero, while the backgrounds extend over a 
wide range. The signal and background PDFs are deter- 
mined from reconstructed MC event samples. The fits are 
shown in Figs. |3Ic) and^d); the signal yields are listed 
in TableQ] We determine B(D~ — v t ) from the e~ and 
H~ samples using Eg. ([3]) and accounting for the decay 



KKtt ' 



fractions of the r 14| . The values obtained are listed in 



Table Q] and are consistent with the previous BaBAR re- 
sult ■ The error- weighted average pj} of the branch- 
ing fractions is B(D~ -> t~v t ) = (5.00 ± 0.35(stai) ± 
0A9(syst)) x 10~ 2 . The weights used in the average 
are computed from the total error matrix and account 
for correlations. As a test of lepton flavor universality 
we determine the ratio B(D~ — > r~ D T ) / B(D~ — » (i~ D^) — 
(8.27 ± 0.77(stat) ± 0.85(syst)), which is consistent with 
the SM value of 9.76. 

As a cross-check of this analysis method, we measure 
the branching fraction for the hadronic decay Dj — > 
K~ K + ir~ . Within the inclusive Dj sample, we re- 
quire exactly three additional charged particle tracks 
that do not overlap with the DKX~f candidate. PID 
requirements are applied to the kaon candidates. The 
mass of the K~K + n~ system must be between 1.93 and 
2.00 GeV/c 2 , and the CM momentum above 3.0 GeV/c. 
We combine the K~K + ir~ system with the signal 7 
and extract the signal yield from the m(K Kir~/) dis- 
tribution. For this mode we choose the loose selection 
m r (DKXj) > 1.82 GeV/c 2 , because this variable is 
correlated with m (K K 717); this corresponds to an in- 
clusive Dj yield of N Ds = (108.9 ± 2.4) x 10 3 . We 
model the signal distribution using reconstructed MC 
events that contain the decay chain D*~ — > D~j and 
D~ — > K~ K + tt~ . In the generic MC and a high statis- 
tics control data sample (for which the inclusive recon- 



Thc weighted efficiency ratio 
is found to be £#-#-^=29. 5%, and we obtain B{D~ — s- 
K-K+iT-) = (5.78±0.20(stat)±0.30(syst))%. The first 
uncertainty accounts for the statistical uncertainties as- 
sociated with the inclusive D~ sample and Nkk-k- The 
second accounts for systematic uncertainties in the signal 
and background models, and the inclusive D~ sample, 
as well as the reconstruction and PID selection of the 
K~K + n~ candidates. This result is consistent with the 
value (5.50 ± 0.23 ± 0.16)% measured by CLEO-c [H[. 

Using the leptonic branching fractions measured 
above, we determine the D~ decay constant using Eq. (fTj) 
and the known values for me, tud,, \V u d\ (we assume 
\Vcs\ = |Kzd|), and the D~ lifetime obtained from 
Ref. [HI]. The fr> 3 values are listed in Table Q] the sys- 
tematic uncertainty includes the uncertainties on these 
parameters (1.9 MeV). Finally, we obtain the error- 
weighted average Jd 3 = (258.6 ± 6A(stat) ± 7.5(syst)) 
MeV. 

In conclusion, we use the full dataset collected by the 
BABAR experiment to measure the branching fractions 
for the leptonic decays of the Dj meson. The measured 
value of Jd b is 1.8 standard deviations larger than the 
theoretical value @, consistent with the measurements 
by Belle and CLEO-c @,[T|. Further work on this sub- 
ject is necessary to validate the theoretical calculations 
or to shed light on possible NP processes. 
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